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Matrix algorithm for determination of the elementary 
paths and elementary circuits using exotic semirings 

Gheorghe IVAN 

Abstract. We propose a new method for determining the elementary paths and elemen- 
tary circuits in a directed graph. Also, the Hamiltonian paths and Hamiltonian circuits are 
enumerated. 

1 Introduction 

Idempotent mathematics is based on replacing the usual arithmetic operations 
with a new set of basic operations, that is on replacing numerical fields by idempo- 
tent semirings. Exotic semirings such as the max-plus algebra R ma z or concatenation 
semiring have been introduced in connection with various fields: graph theory, 

Markov decision processes, language theory, discrete event systems theory, see [3J, [2], 
®- ®- 

In this paper we will have to consider various semirings, and will universally use 
the notation ©, (g>, e, e with a context dependent manning (e.g. := max in R ma3; but 
:= U in ?(£*), e := -co in R max but e := in ?(£*)). 

In many fields of applications, the graphs are used widely for modelling of prac- 
tical problems. This paper will focus on two algebraic path problems, namely: the 
elementary path problem (EPP) and the elementary circuit problem (ECP). 

For a directed graph G = (V,E) (\V\ = n) and u, v G V, (EPP) and (ECP) are 
formulated as follows: 

• (EPP) enumerate the elementary paths from u to v of length k (I < k < n — 1); 

• (ECP) enumerate the elementary circuits starting in u of length k (1 < k < n). 

For to solve the above problems we give a method based on a n x n matrix with 
entries in the semiring of distinguished languages. 

These algebraic path problems are applied into large domains: combinatorial opti- 
mization, traffic control, Internet routing etc. 

The paper is organized as follows. In Section 2 we construct a special idempotent 
semiring denoted by "P*(T,^ W ) and named the semiring of distinguished languages. The 
semiring of matrices with entries in CP*(S^ U) ) is presented in Section 3. This algebraic 
tool is used to establishing of a one-to-one correspondence between the distinguished 
words and the elementary paths in a directed graph. In Section 4 we give an algorithm 
to determine all the elementary paths and elementary circuits in a directed graph. This 
new practical algorithm is based on the latin composition of distinguished languages. 
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Another method for determination of the elementary paths is the well-known latin 
multiplication technique of Kaufmann (see [3] ) . 

2 Semiring of distinguished formal languages 

We start this section by recalling of some necessary backgrounds on semirings 
for our purposes (see j5J, p], [5] and references therein for more details). 

Semirings. Let S be a nonempty set endowed with two binary operations, addi- 
tion (denoted with ffi) and multiplication (denoted with (8)). The algebraic structure 
(S, ©, <8>, e, e) is a semiring, if it fulfills the following conditions: 

(1) (S, ©,e) is a commutative monoid with e as the neutral element for ©; 

(2) (S, <X>, e) is a monoid with e as the identity element for ®; 

(3) <g) distributes over ffi; 

(4) e is an absorbing element for (g>, that is ag)£ = e^a = £, Va£5. 

A semiring where addition is idempotent (that is, a © a = a, Va £ S 1 ) is called an 
idempotent semiring. If & is commutative, we say that 5 is a commutative semiring. 

In the following we introduce the monoid of distinguished words over an alphabet. 

An alphabet is a finite set £ of symbols. An word over £ is a finite sequence of 
symbols of the alphabet £. The number of symbols of a word x is called the length of 
x and its length is represented by \x\. The empty word, denoted with A, is the word 
with length zero (that is, it has no symbols). 

A simple word of length k with 1 < k < n over the alphabet £ (|£| = n) is a word 
of the form a = aia2--.ak~iak such that Oj ^ aj for i ^ j and i,j = 1, k. 

A simple cyclic word of length k + 1 with 1 < k < n over £ is a word of the form 
6 = ai<i2 . . . a^-iai-ai such that a = a±a2 ■ ■ ■ ak-iak is a simple word. The simple cyclic 
words of length 2 over E are the words of the form u = aa for a £ E. 

By a distinguished word over E we mean a word w over E which is a simple word 
of length k or a simple cyclic word of length k + 1 with 1 < k < n. Denote the set of 
distinguished words over E completed with the word A by E^ w . 

For example, if E = {a, 6}, then E^ = {A, a, 6, aa, bb, ab, ba, aba, bab}. 

It is easy to prove that, if E is an alphabet with |E| = n symbols, then E^ is a 
finite set with a n elements, where 



On the set E^ we introduce the binary operation given as follows: 
(1) for all distinguished word x G E^ and simple cyclic word c € E^, we have 



(2) Let x = a\a2 ■ ■ ■ ak-ia^ and y = 6162 . . . o r _i6 r be two simple words of the lengths 
k and r with 1 < k,r < n. The word x o^y £ E^ w is defined by: 




n-l 



(2.1) 



Ao^x = £o^A = A and c x = x c = A; 



(2.2) 
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a 1 ...a k b 2 ...b r ^ib r , if a k = bi, {ai, . . . , a k } (~1 {b 2 , . . . , b r } = 
xo^y = I ai...a k b 2 ...b r -iai, if a k = bi, b r = a±, {ai, . . . , a k } D {b 2 , ■ ■ ■ , b r - 1} = 
A, otherwise. 

(2.3) 

The operation o e is called the latin composition of distinguished words. 

Example 2.1. The set Y? dw of distinguished words over E = {1,2,3,4} has cr 4 = 129 
elements. If x = 123, y = 31, 2 = 1, Ci = 22 and C2 = 343, then 

ci oj 2; = 22 o ; 123 = A, xo l c 2 = 123 343 = A, x o, j/ = 123 o, 31 = 1231, 

t/0(X = 31oi 123 = 3123, z o l x = I 01 123 = 123. Note that xo iy ^yo lX . □ 

For all a, b, c G £^ we have (ao^fe) o;c = aof (bo^c). Then o^, A) is a monoid, 
called the monoid of distinguished words generated by alphabet S. 

A distinguished language L over alphabet S is a subset of distinguished words over 
E completed with the empty word A, that is L C Ti* dw and A € L. 

Convention, (i) If L = {A} or L = {a} where a is a distinguished word, then we 
shall use the notations A and a, respectively. 

(ii) If L 7^ {A}, then we enumerate only its distinguished words of length k > l.D 

The monoid S , contains two special distinguished languages, namely: the language 
{A} (it contains only the empty word A) and the language S (it contains all symbols of 
the alphabet and the empty word A). 

The set of distinguished languages over an alphabet S is ^(E^). 

Since distinguished languages are sets, all the set operations can be applied to dis- 
tinguished languages. Then the union and intersection of two distinguished languages 
are distinguished languages. 

The latin composition of the distinguished languages L\ and L 2 is the distinguished 
language L\ o e L 2 defined by 

Li o?L 2 = {x o £ y I x € Li and y £ L 2 }, (2.4) 

that is, L\ 0£ L 2 is the set of distinguished words obtaining by the latin composition of 
words of L\ with those of L 2 . 

For example, if L\,L 2 C E^ w are distinguished languages over E = {1,2,3,4}, 
where L x = {2,412} and L 2 = {11,23}, then 

L x o e L 2 = {2, 412} o e {11, 23} = {2 o t 11, 412 o t 11, 2 o e 23, 412 o e 23} = {23, 4123}. 
Proposition 2.1. (i) Let L, L 1: L 2 , L 3 e ^(E^J. Then: 
L 0l A = A 01 L = A, Lo,S = SojL = L, (Li o z L 2 ) o ; L 3 = L x (L 2 L 3 ). (2.5) 
(ii) T/ie set "P*(T, dw ) endowed with multiplication 01 has a structure of monoid. 
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Proof, (i) Using the definitions one easily verify that (2.5) holds. 

(ii) From (i) it follows that (P* (Y,* dw ) , o h £) is a monoid. □ 

On the set T*(E^) of distinguished languages we define the binary operations: 

Li0L 2 :=LiUL 2 and L x © L 2 := L x o e L 2 , VL X , L 2 € 3>*(E2J. (2.6) 

Proposition 2.2. (E dw ) , U, o^, A, S) is an idempotent semiring. 

Proof. For to verify the conditions from definition of an idempotent semiring we 
apply the properties of the union of sets and Proposition 2.1. □ 
We call (IP*(E^ W ), U, o^, A, X) the semiring of distinguished languages over S. 



3 Matrices over semirings and directed graphs 

Let (S, ©,©,£, e) be an (idempotent) semiring. For each positive integer n, let 
M n (S) be denote the set ofnxn matrices with entries in S. The operations © and 
© on S induce corresponding operations on M n (S) in the obvious way. Indeed, if 
A = (Aij),B = (Bij) G M n (S) then we have: 

A(B B = ({A © B)ij) and ^ © B = ({A © i,j = T~n where 

n 

(A © B)ij := ^ © B - and (A © := A ik © J3jy . (3.1) 

k=i 

The set M n (S) contains two special matrices with entries in S, namely the zero 
matrix 0©„, which has all its entries equal to e, and the identity matrix I® n , which 
has the diagonal entries equal to e and the other entries equal to e. 

It is easy to check that the following proposition holds. 

Proposition 3.1. (M n (S), ©, ©, 0® n , I® n ) is & n idempotent semiring, where the op- 
erations © and © are given in (3.1). □ 

We call (M n (S), ©, ©, 0© n , I®n) the semiring ofnxn matrices with entries in S. In 
particular, if S := (IP*(£*J, U, o e , 0, £), then (M n (?*(£*J), U, o*, O e „, J® n ) is called 
the semiring of n x n matrices over 3>*(Ej ). The operation © := is called the 
multiplication of matrices based on latin composition of words. 

Example 3.1. Let be semiring (M 2 (T (E* dw )) , U, o e , O e2 , 1®2) with £ = {a,fe,c}. The 
product i Of B of the A, B with entries in the semiring J > *(S^ U) ) is 

e \ bca be J \ c e J \ (bca og b) © (fee og c) (feca og ab) © (fee e) 

afeffie effie \ _ / a&Ue eUe A _ a & e i 
£ © fee fecafe © e) ) ~ I e U fee fecafe U e) / ~~ I fee fecafe 
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A directed graph is a pair G = (V, E) where V is a finite set of vertices of the graph 
G and E Q V x V is & set of arcs of G. A typical arc (u, v) G E is thought of as an 
arrow directed from u to v. 

Let G = (V, E) be a directed graph with \V\ = n. A pai/i /rom u to v of length k 
(k > 1) in G is a sequence of vertices p = (ui,i>2, • • • ,^,^+1) with v\ = u, Vk+i = v 
such that € E for all i = l,k; v\ is called the starting vertex and Vfc+i the 

end-vertex of p, respectively. The length of path p will denoted by £{p). 

A path p = (vi, V2, ■ ■ ■ , Vk, ffc+i) is called circuit if v^+i = t>i and fc > 1. In 
particular, for A; = 1 we obtain the circuit (t>i, v\) of length 1. 

We denote with P(vi,Vj, k) (k > 1), the set of all paths of length k from the starting 
vertex V{ G V to end-vertex «j G V. In particular, when Vi = Vj, C(vi, k) = P(vi, Vi, k) 
( k > 1) is the set of all circuits of length k starting at vertex Vi. 

A path p = (vi, i>2, • • • , Vk, Vk+i) is called an elementary path from v\ to Vk+i, if 
A; > 1 and Vi ^ for j ^ j and i,j = 1, fe + 1. A circuit c = (ui, U2, • • • , ttfc, tti) with 
£(c) = k is called an elementary circuit, if («i,U2, . . . , u^) is an elementary path. 

We denote with P e i em (vi,Vj,k) (k > 1), the set of all elementary paths of length k 
from Vi £ V to Vj G V. In particular, when = Vj, then C e i em (vi, k) = P e iem{vi, Vi, k) 
(k > 1) is the set of all elementary circuits of length k starting at Uj. 

A Hamiltonian path (resp., circuit) is a path (resp., circuit) that contains each 
vertex exactly once. Hence, a Hamiltonian path (resp., circuit) is an elementary path 
Ph with £(ph) = n — 1 (resp., an elementary circuit ch with £(ch) = n). 

A weighted directed graph is a graph G = (V, E) with a mapping w : E ^ S that 
assigns each arc (u, u) G E a weight w(u, u) from the semiring (S 1 , ©, (g), e, e). A weighted 
directed graph with the cost function w is denoted by G = (V, £7, u;). 

The weight or cost of path p = (v±, . . . , Vk,Vk+i) is the element w(p) G S where 



To each given weighted directed graph G = (V, E, w) with V = {vi,V2, ■ ■ ■ , v n } we 
can associate a n x n matrix M W (G) with entries in a semiring (S, ©, 0, e, e) as follows. 
For this, we define the matrix M W (G) = (My) G M(n, S) where 



To each directed graph G = (V, E) with V = {v\,v%, . . . , v n } we can associate two 
weight functions in the following way. 

• Let be the numbering semiring (N, +, •, 0, 1) of natural numbers endowed with 
the usual addition and multiplication. Consider the weight function w a : E — >■ N 
defined by w a (vi,Vj) = 1 for all (yi,Vj) G E. The matrix M Wa (G) G M(n,N), denoted 
with A, is called the adjacency matrix of graph G. 



k 



w (p) = <§Qw(vi,v i+1 ). 



(3.2) 



i=l 




(3.3) 
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• Let be the idempotent semiring (J > *(E^ UI ), U, o^, 0, S) of distinguished languages 
over alphabet £ = V. Define the weight function we : E — >■ ^(E^) given by 
Wf(vi,Vj) = ViVj for all (vi,Vj) G £7 (that is, we(vi,Vj) is the distinguished lan- 
guage which contains only the distinguished word fjfj of length 1). The matrix 
M We (G) € M(n, J , *(SJ )) is denoted with L and is called the Zaim matrix of G. 

More precisely, the adjacency matrix A = (Aij) 6 M(n, N) and latin matrix L = 
(Ljj) € M(n, y*(S^)) associated to graph G are given by 

^ = (! * -d ^=(^ * (3.4) 

J [ if {Vi,Vj)£E { £ lf (Vi,Vj)pE 

4 Matrix algorithm for enumerating of elementary paths 
and elementary circuits in a directed graph 

Consider the latin matrix L = (Lij) € M(n, y*(S^)) associated to graph G = 
(V,E) (\V\ = n), defined by (3.4). Using the multiplication of matrices based on latin 
composition, we define by recurrence the power of matrix L in the following way: 

LW=L, L [2] =Lo e L, L [k] = L o e L [k - 1] for k > 2. (4.1) 

Applying (3.1) and replacing © and (g> with the correspondent operations of the 
semiring M(n, y*(S^ w )), we have = (L^) for i,j = l,n, where 

n 
m=l 

It is easy to prove that: 

L [n] is a diagonal matrix and L [n+q] = 0© n for all q > 1. (4.3) 

Theorem 4.1. Lei G = (V,i2) 6e a directed graph with V = {v\, V2, ■ ■ ■ , v n } and the 

latin matrix L = (Lij) € M(n, T*(Z* d J). If = (L [ ^), i/ien: 

4j = P e lem(vi,Vj, k), i,j = T~n, i ^ j, 1 < k < n - 1; (4.4) 

^li' = Ceiem(fi, fc), Z = 1, n, 1 < k < U. (4.5) 

Proof. We proceed by induction on k. Choose Vi and Vj arbitrarily. 

(1) Case i ^ j. For k = 1, the relation (4.4) holds. Indeed, we have that the 
only elementary path in P e i em (vi,Vj, 1) is the arc (vi,Vj), since = VjWj. Therefore, 
Pelem(vi,Vj,l) = {(vi,Vj)}. Note that if = e, then P e i em (vi,Vj,l) = 0. 

Now assume the relation (4.4) holds true for k. Then L^j = P e iem(vi, Vj, k), that is 
llfj represent the set of all elementary paths of length k from Vi to Vj . 
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Using (4.2) the entry of the matrix =Lo f is given explicitly by 

= (L il o e L^)U...U(L im o e L^ ] j )U...U(L in o e L^ ] j ), k > 1. (4.6) 

We first evaluate the term L im o £ Z,j^l from the equality (4.6), for an fixed integer 
m with 1 < m < n. We have the following situations. 

• If Li m = e or L^jj = e, then Lj m o £ L m ^ = e, that is the set of elementary paths 
of length k + 1 from Vi to stopping at vertex v m is empty set. 

r[k] 

J mj 



• Li m 7^ ^ and iM- / e. By the induction hypothesis we have 



L mj = p eiem{v m , Vj,k) = {p m jiP m j, ■ ■ ■ ^m/'Pmj}' where pj^. is an elementary path of 
length k from v m to Vj for 1 < s < r. Since L« m = ViV m it follows that 

^im °i ifi = ViV m Oi {p mj ,p mj , ■ ■ ■ ,Pmj~ iPmj) 

where p m j is regarded as a distinguished word, having the set of vertices 

x Lj = {"» = v mod ,v mo+ ij, . . . ,v mo+ k~i,j,v mo+ k,j = vj}. The element ViV m °e P mj 
can be take the following values: 

- ViV m 0£ p S m j = {ViV m Vmo+l,j ■ ■ ■ Vm +k-l ,jVj }, if {v^} (~l X mj = 0, 

that is ViV m °£P m j ^ s an elementary path of length k + 1 from to u^; 

- Ui«m Pmj = if {«»} n ^mj + 0- 

Then Lj m L^- = {A} or Lj m L^- is a set which contains t (1 < t < s) 
elementary paths of length k + 1 from to stopping at the vertex v m . 

Therefore, L [ ^ +1] = (J (L, m o £ lJJ) = P e i em (vi,Vj,k + 1). Hence, (4.4) holds for 

m— 1 

k + 1. This completes the inductive step and proves the assertion in the case i ^ j. 

(2) Case i = j. If = vm, then C e z em (uj, 1) = {(vi,Vi)}. Also, if I,** = e, then 
C e i em (vi, 1) = 0. Hence, (4.5) holds for k = 1. Assume that (4.5) holds true for fc. Then 

[fel / \ [fel 

^ii = C e i em (Vi,k), that is L^ J represent the set of all elementary circuits of length k 
starting at Vi. 

The entry of the matrix L^ is lJJ +1] = Q (L im o e Z,jJJ), k > 1. 

m=l 

We evaluate the term Lj m for an fixed integer m with 1 < m < n. We have 

the following situations. 

[fel [fel 

• If L irn = e or L m \ = e, then L im °e L m \ = e. 

• L im + £ and L m \ ^ e. Applying (4.4) we have L^ ] - = P elem {v m ,Vi, k) = 
illnv Imv •••■> Im^jQmi}, where q mi is an elementary path of length k from v m to Vj for 
1 < s < r. Since Li rn = ViV rn it follows that 

Lim °l L l Si = ViV rn °£ {q m i, Qmi, Qm^Qmi) 

where q mi is regarded as a distinguished word, having the set of vertices 
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Y mi = i v m = v mu i,v mi+lji , . . . , v mi+k -i t i, v mi+kii = Vi}. The element viv m o e q s mi can 
be take the following values: 

- ViV m o e q s mi = {viV m v mi +i,i ■ ■ ■ u mi+ jfe_ 1( i«i}, if {v^ n (Y^ \ {vi}) = 0, 

that is ViV m o£ q s mi is an elementary circuit from Vi to Vi of length k + 1; 

- v lVm o e q* mi = {A}, if { Vi } n (F^ \ {«<}) ^ 0. 

Then Lj m L^l = {A} or Lj m is a set which contains t% (1 < h < s) 

elementary circuits of length k + 1 from to V{ stopping at v m . Therefore, L^ +1 ^ = 
C e i em (vi,k + 1). Hence, (4.5) holds for k + 1. This completes the inductive step and 
proves the assertion in the case i = j. □ 

Applying Theorem 4.1 we give an answer of (EPP) and (ECP) for a directed graph. 
In this purpose we give a new method based on latin composition of distinguished 
languages, We will called the algorithm of latin composition of distinguished languages 
(shortly, LCDL-algorithm). 

For a directed graph G = (V,E) with V = {v\,V2, ■ ■ ■ ,v n }, the LCDL- algorithm 
consists from the following steps: 

Step 1. Associate the latin matrix L E M(n,3 5 *(S^)) to graph G; 

Step 2. For each k(l < k <n) compute the matrix = (L^); 

Step 3. (i) For each pair and each k(l < k < n — 1) enumerate the elementary 
paths of length k from V{ to Vj in G (we apply the relation (4.4)); 

(ii) For each i and each k(l < k < n) enumerate the elementary circuits of length 
k starting at Vi in G (we apply the relation (4.5)). 

Remark 4.1. Let G = (V, E) be a directed graph with V = {vi,V2, ■ ■ ■ , v n }. 

(1) (1) For 1 < k < n — 2 and i,j = l,n with i ^ j, the elements indicates the 

elementary paths which are formed by k arcs, so that l|™ ^ determines all Hamiltonian 
paths in G between u$ and Vj. 

(2) For 1 < k < n — 1 and i = l,n, the elements indicates the elementary 

In] 

circuits of length k, so that L u determines all Hamiltonian circuits in G starting at V{. 

(ii) The necessity to determine the elementary paths and elementary circuits of 
maximum length in a directed graph arises. The LCDL-algorithm determines the 
Hamiltonian paths (resp., Hamiltonian circuits) when there exist or determines the 
elementary paths (resp., circuits) of maximum length when we have no Hamiltonian 
paths (resp., circuits). □ 

Remark 4.2. (Jiyj The powers of the adjacency matrix A £ M(n; N) associated to 
graph G = (V, E) are used to find the number of distinct paths and distinct circuits 
between two vertices in V (two paths or circuits of length k are distinct if they visit a 
different sequence of vertices). More precisely: 

Let A k = (Afj) be the k—th power of the adjacency matrix A. Then: 

= \P(vi,Vj,k)\ and A^ = \C(v h k)\ for all k>l. □ 
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Example 4.1. Let G = (V,E) be a directed graph with vertex set V = {vi, v 2 , v%, v 4 } 
and the adjacency matrix A G M 4 (N) of G where 

/ 1 1 1 1 \ 

A= 0111 

1 ' 

V / 

The geometric representation of graph G is given in Figure 4.1. 



Q 




v3 



Fig.4.1. 



Computing the power A 3 of the adjacency matrix A in terms of numbering semiring 
(N, +, ■, 0, 1) we have Af A = 5 and A\ 2 = 1. Then: 

• Between v\ and V4 there exist 5 paths of length 3, namely: p\ = (vi,V2,V3,V4), 



Pi = (vi,t>l,Vl,V4), Pi = (V1,V 1 ,V 2 ,V4). 



P\ 



{V1,V 2 ,V 2 ,V4) 



pI 



(vi,wi,U3,w 4 ); 

• G has one circuit of length 3 starting at v 2 , namely cf = (^2,^2,^2,^2)- 
(ii) We consider the alphabet S = {v±, v 2 , V3, v^}. For to find the elementary paths 
and elementary circuits in graph G we apply LCDL-algorithm. 
The latin matrix L € M^J^E^)) associated to graph G is 



L 



e 

£ 
£ 



V\V 2 
V 2 V 2 

£ 

£ 



V1V3 
V 2 V 3 

£ 

£ 



V1V4 \ 

V 2 V 4 

V3V4 

Z ) 



We compute the powers od the latin matrix L in terms of the semiring of distinguished 
languages. We first compute the matrix = We have 



L [2] = L Of L 



( £ £ ViV 2 V 3 {V!V 2 V4,V1V 3 V4} \ 

£ £ £ V 2 V3V4 

£ £ £ £ 

\ £ £ £ £ J 
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For example, L\l is computed as follows 
L\l = (L u o e L i4 ) U (L12 L 24 ) U (L13 L 34 ) U (L14 o £ L 44 ) = (i^i o e u 1 u 4 )U 

U(uiW2°^W2V4)U(wif30£ti3ti 4 )U(t'iW40^e) = £U{f iV2f4}U{wiW3W 4 }Ue = {f 1V2V4, V1V3V4}. 
[21 

Since L^ 4 J = {viV2^4, V1V3V4} it follows that there exist two elementary paths of 
length 2 from v\ to w 4 and we have P e iem(v 1, t>4, 2) = {(^1,^2,^4), (^1,^3,^4)}. 
The matrices i' 3 ' = L o f il 2 l and = L op L^ are given by 



L [3] 



V 





£ 


£ 


V\V2 ^3^4 


£ 


£ 


£ 


£ 


£ 


£ 


£ 


£ 


£ 


£ 


£ 


£ 



and LM 



f £ S S £ \ 

£ £ £ £ 

£ £ £ £ 

\ £ £ £ £ / 



Using the matrices and one obtains the following results: 

[31 

• L^ 4 J = {viV2^3f4}- Then P e j em (ui, U4, 3) = {(f 1, t>2, v 3 , v 4 )} and G has only one 
Hamiltonian path. The set of elementary paths of maximum length from V2 to V4 is 

[21 [31 

P e iem{v2,V4,2) = {(v 2 , v 3 , V4)}, since L 24 = {V2V3V4} and L\l = e. 

• l]^ = £ for 2 < fc < 4 and i = 1,4. Then C e i em (vi,Vi, k) = 0. The elementary 
circuits of maximum length are those of length 1. □ 



Application. The finding of Hamiltonian paths and Hamiltonian circuits 
of minimal cost in a weighted directed graph. 

The LCDL- algorithm can be to use for solving of the following two problems in a 
weighted directed graph G = (V, E) : 

(i) find a Hamiltonian path of minimal cost between two vertices in G; 

(ii) for each v £ G, find a Hamiltonian circuit of minimal cost starting at v. 

One way to solve the above problems consists of searching all possible Hamiltonian 
paths and Hamiltonian circuits (we apply LCDL-algorithm) and computing their cost 
(we apply the relation (3.1)). 



Example 4.2. Let G = (V, E, w) be a weighted directed graph where V = {1,2,3, 4, 5}, 

E = {(1, 2), (1, 3), (1, 5), (2, 1), (2, 5), (3, 2), (4, 3), (4, 5), (5, 1), (5, 2), (5, 3), (5,4)} and the 
cost function w cost : E -> R, (i,j) >-> w cost (i,j) = Wij given by 



(i,3) 


(1,2) 


(1,3) 


(1,5) 


(2,1) 


(2,5) 


(3,2) 


(4,3) 


(4,5) 


(5,1) 


(5,2) 


Wij 


4 


2 


6 


3 


3 


1 


5 


4 


6 


1 


(iJ) 


(5,3) 


(5,4) 




Wij 


2 


1 



Matrix algorithm for determination of the elementary paths and 
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The latin matrix L G M 5 C?(T,* dw )) (S = {1, 2, 3, 4, 5}) associated to graph G is 



/ e 12 13 

21 £ £ 

£ 32 £ 

£ £ 43 



15 \ 

25 

£ 

45 



V 51 52 53 54 e J 
Let us we compute the powers of the latin matrix L. We have 



L [2] = L oc L 



( {121,151} {132,152} 153 154 

251 {212,252} {213,253} 254 

321 £ £ £ 

451 {432,452} 453 454 

\ 521 {512,532} {513,543} e 



125 
215 
325 

£ 

{515,525,545} / 



The matrix = L o e L^ 2 ' has the following lines: 



Ll 



[3] 



L 



({1321,1521,1251}, 1532, {1543,1253}, 1254, 1325); 

(e, {2512,2132,2532,2152}, {2513,2543,2153}, 2154, e); 

(3251, £, {3213,3253}, 3254, 3215); 

({4321,4521}, {4512,4532}, 4513, e, 4325); 

(5321, {5132,5432}, 5213, e, {5215,5125,5325}). 



The matrix = L has the following lines: 



L 



L 



L 



L 



L 



({15321,13251}, 15432, 12543, 13254, e); 
(£, {25132,25432,21532}, 21543, e, e); 
(£, £, {32513,32543,32153}, e, e); 
({45321,43251}, 45132, 45213, 43254, 43215); 
(54321, £, £, £, {53215,51325,54325}). 



Finally, we have 



L( 5 ) = o £ L 



( 154321 

£ 
£ 
£ 

£ 



£ 

215432 

£ 
£ 
£ 



£ 
£ 
£ 

432154 

£ 



543215 / 

• From it follows that G has 11 Hamiltonian paths. For example, P e i em (4, 1, 4) = 
{p\ H = (4,5,3,2,1),^^ = (4,3,2,5,1)} since lJJ = {45321,43251}. We have 
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w (pf h) = 10 and w{p2 H ) = 15. It follows that p\ H is a Hamiltonian path between 
the vertices 4 and 1 having the maximal cost equal to 15. 

• From it follows that G has 4 Hamiltonian circuits. For example, C e i em (l, 1, 5) = 
l[ 5 J = {cf ^ = (1,5,4,3,2,1)}. We have u>(cf H ) = 16. Hence c\ H is a Hamiltonian 
circuit starting at vertex 1 having the maximal cost equal to 16. □ 

Conclusions. The LCDL- algorithm can be easily programmed and gives an 
efficient solution of EPP and ECP for a finite directed graph. This can be seen as an 
improved version of Kaufmann's algorithm. 

Let us list some classes of algebraic path problems, which can be reduced to applying 
of the LCDL-algorithm: 

(i) enumeration of the elementary paths (resp., circuits); 

(ii) determination of the elementary paths (resp., circuits) of maximum length; 
(Hi) testing a graph G for having Hamiltonian paths or Hamiltonian circuits; 

(iv) optimization (Hamiltonian path or Hamiltonian circuit of minimal cost). □ 
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